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Abstract 

The event-by-event distributions of harmonic flow coefficients v„ for n-2-A are measured in Pb-Pb collisions at 
y*AW = 2.76 TeV with the ATLAS detector, using charged particles with pj > 0.5 GeV and \rj\ < 2.5. The shape of 
the v„ distributions is consistent with Gaussian fluctuations in central collisions for and over the measured centrality 
range for V3 and V4. When these distributions are rescaled to have the same (v„), the resulting shapes are nearly the 
same for pj > 1 GeV and 0.5 < pj < 1 GeV. The v„ distributions are compared with the eccentricity distributions 
from two initial geometry models: Glauber and the MC-KLN model. Both fail to describe the V2 data consistently 
over the full centrality range. 



In recent years, the measurement of harmonic flow coefficients v„ has provided important insight into the hot and 
dense matter created in heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron 
Collider (LHC). These coefficients are generally obtained from a Fourier expansion of particle azimuthal angle distri- 
butions, ^ oc 1 + 2 2™ =1 v„cosn(0 - <£„), where 0„ represents the phase of v n (event plane or EP) Q]|2). Previous 
measurements determined the v n from the distribution of - 0„, accumulated over many events. This event-averaged 
v„ mainly reflects the hydrodynamic response of the created matter to the average collision geometry in the initial 
state. However, more information can be obtained by measuring the v„ on a event-by-event (EbE) basis. This is par- 
ticularly useful for understanding the nature of the EbE fluctuations in the initial collision geometry. This proceedings 
present the methods and results of EbE v„ for n = 2 — 4 obtained with the ATLAS detector (3). 

This analysis is based on 8 pb~ l of minimum bias Pb-Pb data collected in 2010 at -Js^ - 2.76 TeV [4|. The v„ 
coefficients are calculated using tracks with pj > 0.5 GeV and \rj\ < 2.5, reconstructed by the inner detector (ID). 
To illustrate the level of EbE fluctuations, Figure [T] shows the distributions of track <p and track pair relative angle Atp 
with pj > 0.5 GeV for three central events. Rich EbE patterns, beyond the structures of the detector acceptance (solid 
points) and statistical fluctutions, are observed. These distributions are the inputs for the EbE v„ analyses. 

Two methods are used to obtain the EbE v„. The first method starts with a Fourier expansion of the azimuthal 
distribution of charged particles in a given event, reweighted by the inverse of the tracking efficiency for each particle: 



— - 1 + 2 £ vf 8 cos n0 - ®° bs ) = 1 + 2 £ (v* cos wj> + v£ sin #)) , v° ta = yj (i#) 2 + (<f , tann*f = (1) 

where v° bs is the magnitude of the observed per-particle flow vector: v w obs = (v° b x s , v° b y s ). In the limit of infinite 
multiplicity and the absence of non-flow effects, it approaches the true flow signal: v° bs — > v„. The key of the 
measurement is to determine the response function p(v n obs |v„) or p(v° bs |v„), used to unfold these smearing effects. 

In order to determine the response function, the tracks in the ID are divided into two subevents with symmetric 77 
range, rj > and 77 < 0. The smearing effects are estimated from the difference of the flow vectors between the two 
subevents, /? su b ((v„ obs ) a - (v„ obs )''), for which the physical flow signal cancels. This distribution is observed to be well 
described by a 2-D Gaussian with identical widths in both dimensions, 6 2SE , and hence the 2-D response function can 
be obtained as p(v° bs \v n ) = cp suh (c [(v„ ob T - (v„ obs )'']), with c = 2 and V2 for the full-ID and half-ID, respectively. 



'A list of members of the ATLAS Collaboration and acknowledgements can be found at the end of this issue. 
Preprint submitted to Nuclear Physics A 



September 20, 2012 




40000 


" Run 169927 


ATLAS Preliminary 




Event 400484 


Pb-Pb^=2.76TeV 




0-5% | | 


P t >0.5 GeV 


39000 


















£ 38000 






37000 








Figure 1: Single track 
<p (top) and track pair A<f> 
(bottom) distributions for 
three events (from left to 
right) in 0-5% centrality in- 
terval (4). The pair dis- 
tributions have been folded 
into [-0.5;r,1.5/r]. The bars 
indicate the statistical un- 
certainties, the solid curves 
indicate a Fourier param- 
eterization including first 
six harmonics: dN/d(f> = 
A(l + 2 2* =1 c„cosn((p - 
¥„)) f° r single track dis- 
tributions and dN/dA<p = 
A(l + 2 2f =1 c„cosn(A0)) 
for pair distributions, and 
the red solid points indi- 
cate the acceptance func- 
tions (arbitrary normaliza- 
tion). 



A simple shift of this 2-D distribution to v„ = (v„ x , v„ y ), followed by an integration over the azimuthal angle, gives 
the desired 1-D response function in the form of the well-known Bessel-Gaussian function 0: 



P(v° b >„) . 



v° b V 



"/», 



6 2 



Jy/2 for half-ID 
J2 for full-ID 



(2) 



where Iq is the modified Bessel function of the first kind and the difference s = v° - v„ accounts for the smearing. 

In the second method, the observed flow signal is defined from an EbE pair distribution, obtained by convolving 
the tracks in the first half-ID with those in the second half-ID: 



dN 
dAcp 



where S\ 



1+2 ^ (v° b x Sl cos n(pi + vjfy 1 sin n(pi ) ® 1+2 ^ (v°J 2 cos nfc + v* 5,2 sin rup^j =1+2 ^ (A„ cos nAcp + B„ sin nA(f>) 

n J L n in 

= (Al + Bl)" 4 = [(v* 1 ) 2 + (vC 1 ) 2 ]" 4 [(v° b ") 2 + (v* 2 ) 2 ]" 4 = V^ bS ' v « bS2 = V(v„ + ^)(v„ + . S2 ) (3) 

v„ are independent variables described by the probability distribution in Eq. [2] 



obsi 



v„ and 52 



with 6 = c\ SE / V2. The response function for v^ bs ' 2PC is different from v° bs due to the presence of two random variables 
The Bayesian unfolding procedure from J6| is used to calculate the v„ distribution, in three p T ranges: p T > 0.5 
GeV, I > pj > 0.5 GeV and pj > 1 GeV. In each case, the prior (initial distribution) is taken as the v° bs distribution 
from the full-ID, and the number of iterations A^i ter is adjusted according to the sample statistics and binning. The 
convergence is generally reached for Ni tel > 8 in the case of n = 2, but more iterations are required for larger n 
and in more peripheral collisions. The final results are found to be independent of the choice of priors, although the 
convergence is faster with an appropriate choice. The robustness of the unfolding is checked by comparing the results 
measured independently for the half-ID and the full-ID, as well as by comparing between the two unfolding methods. 
Despite the large differences in their initial distributions, the final unfolded results always agree. The final results 
are also found to be insensitive to the explicit requirement of a 77 gap between the two subevents in both methods, 
suggesting that the influences of the short-range non-flow correlations on the final results are not significant. 

Figure [2] shows the probability distribution of the EbE v„ in several centrality intervals obtained for charged 
particles at pj > 0.5 GeV. The shape of these distributions changes strongly with centrality for V2, while it is relatively 
unchanged for higher-order harmonics. These distributions are compared with the PDF obtained by radial projection 

of a 2-D Gaussian distribution in v„: P(v„) = ^-e ^ ,cr = y| (v„). The Gaussian description works well for V3 and 
V4 over the measured centrality range, but fails for V2 beyond the top 2% most central collisions. 
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Figure 2: The proba- 
bility distribution of the 
EbE v„ for n = 2 (left), 
« = 3 (middle) and n = 
4 (right) H. The solid 
curves are Gaussian dis- 
tributions with mean ad- 
justed to the measured 
(v„), shown for the 0-1% 
centrality interval for V2, 
but for all centrality inter- 
vals for v»3 and V4, 



Many quantities can be calculated directly from these distributions, such as the mean (v„), width <x,, n , ratio cr Vn / (v n ) 
and RMS value -\f(y^} s ^j{v„} 2 + ~o% n - The cr v J(v n } is a measure of the relative fluctuations of v„ and was previously 
estimated indirectly from the two- and four-particle cumulant methods [7|. Figure [3] shows the cr v J{v n ) calculated for 
the three pj ranges. Despite the strong pj dependence of (v„) and cr v , the ratios are remarkably stable, suggesting 
that the hydrodynamic response to the initial geometry is nearly independent of pj. For V2, the values of tr Vn /(v n ) 
vary strongly with (N pmt ), and reach a minimum of about 0.34 at {N pmt ) ~ 200 or 20-30% centrality range. For V3 and 
V4, the values of cr Vl J{v n ) are almost independent of (N p - drt ), and are consistent with the value expected from Gaussian 
distributions ( ~\J4/tt - 1 * 0.523 as indicated by the dotted lines). 
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Figure 3: The tr v „ /<v„) vs. (N pM ) 
in three pj ranges for n = 2 (top- 
left), n = 3 (bottom-left) and n = 4 
(bottom-right) (4). The dotted lines 
indicate - 1 a 0.523 expected 

for Gaussian fluctuations. Top-right 
panel shows the <x f , /(ft) for the 
Glauber model |8j and the MC-KLN 
model 0. 
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Figure [4] compares the EbE v'2 distributions with the distributions of the eccentricity 62 of the initial geometry, 
calculated for the Glauber model and the MC-KLN model (version 3.46) [9|. The €2 distribution for each centrality 
interval is rescaled to match the (v'2) of the data, and then normalized into a PDF. Figure [4] shows that the rescaled 
62 distributions describe the data well for the most central collisions, but start to fail in non-central collisions. This 
behavior is also reflected in the comparison of <x,, 2 / (V2) with cr ft /{ei) in the top-right panel of Figure|3] The agreement 
with the models for n — 3 - 4 (see @ for more details) are better than the n — 2 case, however, this could simply 
reflect the fact that all distributions are dominated by Gaussian fluctuations, which have a universal shape. 

The EP method in general is known to measure a v„ value between the simple average and the RMS of the true 
v„ IfTUl : (v„) < v^ p < y/{v%). This relation is checked explicitly in Figure [5] based on the EbE v„ distributions. For 
V3 and V4, the values of v„ F are almost identical to ^(v%); For v>2, the values of vf p are in between (v„) and y/{v%): 
they are closer to (v„) in mid-central collisions where the EP resolution factor is close to one, and approach y(v^) in 
peripheral collisions where the resolution factor is small. 
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Figure 4: The EbE vj dis- 
tributions compared with the ei 
distributions from the Glauber 
model (red lines) and the MC- 
KLN model (blue lines) (4). 



Figure 5: Top panels: Com- 
parison of (v n ) and J{\%) = 

^(v„> 2 + o-i with vf |4 1. 
Bottom panels: the ratios of 

A /<^> and vf to <v„> |4|. 
The dotted lines in bot- 
tom panels at J(\%)/(v„) = 

V4/7T a 1.13 indicate the 
value expected for Gaussian 
fluctuations. 



In summary, the EbE v„ distributions for n — 2 - 4 are measured for Pb-Pb collision at yfs^ — 2.76 TeV. The 
shape of the v„ distributions is consistent with Gaussian fluctuation in central collisions (0-2% centrality range) for V2 
and over the full centrality range for V3 and V4. The ratio of the RMS to the mean, cr v J{v n ), is studied as a function of 
(A^ part ) and p T . The values of cr v J(v n ) are found to be independent of p T , suggesting that the hydrodynamic response 
to the eccentricity of the initial geometry has little pj dependence, however they are found to reach a minimum of 
0.34 for v'2 around {N pmt ) ~ 200. A comparison of the v„ distributions with the eccentricity distributions of the initial 
geometry from the Glauber and MC-KLN models, shows that both models fail to describe the data across the full 

centrality range. The v„ from the event plane method is found to lie in between (v„) and ^j{v n } 2 + cr^. These results 
may shed light on the nature of the fluctuations of the created matter in the initial state as well as the subsequent 
hydrodynamic evolution. 
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